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Abstract 

Quantum chromodynamics (QCD) at sufficiently high density is expected to 
undergo a chiral phase transition. Understanding such a transition is of partic- 
ular importance for neutron star or quark star physics. In Lagrangian SU(3) 
lattice gauge theory, the standard approach breaks down at large chemical 
potential //, due to the complex action problem. The Hamiltonian formula- 
tion of lattice QCD doesn't encounter such a problem. In a previous work, we 
developed a Hamiltonian approach at finite chemical potential /i and obtained 
reasonable results in the strong coupling regime. In this paper, we extend the 
previous work to Wilson fermions. We study the chiral behavior and calculate 
the vacuum energy, chiral condensate and quark number density, as well as 
the masses of light hadrons. There is a first order chiral phase transition at 
zero temperature. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is the fundamental theory of strong interactions. It 
is a SU(3) gauge theory of quarks and gluons. Precise determination of the QCD phase 
diagram on temperature T and chemical potential ji plane will provide valuable information 
for the experimental search for quark-gluon plasma (QGP). The ultimate goal of machines 
like the Relativistic Heavy Ion Collider (RHIC) at BNL and the Large Hadron Collider 
(LHC) at CERN is to create the QGP phase, and replay the birth and evolution of the 
Universe. Such a new state of matter may also exist in the core of neutron stars or quark 
stars at low temperature T and large chemical potential ji. Lattice gauge theory (LGT), 
proposed by Wilson [1] is a first principle non-perturbative method for QCD. Although 
it is the most reliable technique for investigating phase transitions in QCD, it is not free 
of problems: complex action at finite chemical potential and species doubling with naive 
fermions. 

In Lagrangian formulation of LGT at finite chemical potential, the success is limited 
to SU(2) gauge theory [2,3], while in the physical SU(3) case, complex action [4,5] spoils 
numerical simulations with importance sampling. Even though much effort [6-8] has recently 
been made for SU(3) LGT, and some very interesting information on the phase diagram at 
large T and small n has been obtained, it is still extremely difficult to do simulations at large 
chemical potential. QCD at large /i is of particular importance for neutron star or quark star 
physics. Hamiltonian formulation of LGT doesn't encounter the notorious "complex action 
problem". Recently, we proposed a Hamiltonian approach to LGT with naive fermions at 
finite chemical potential [9,10] and solve it in the strong coupling regime. We predicted that 
at zero temperature, there is a first order chiral phase transition at critical chemical potential 

= m { ^ n = A4 0) /3, with 

m ^dyn an d M$ being the dynamical mass of quark and nucleon 
mass at /i — respectively. (We expect this is also true for Kogut-Susskind fermions.) 
By solving the gap and Bethe-Salpeter equations, the authors of Ref. [11] obtained the 
critical point the same as ours; but they concluded that the chiral transition is of second 
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order, different from ours. Our order of transition is consistent with other lattice simulation 
results [12]. 

Wilson's approach to lattice fermions [1] has been extensively used in hadron spectrum 
calculations as well as in QCD at finite temperature. It avoids the species doubling and 
preserves the flavor symmetry, but it explicitly breaks the chiral symmetry [1,13-18], one of 
the most important symmetries of the original theory. Non-perturbative fine-tuning of the 
bare fermion mass has to be done, in order to define the chiral limit [19,20]. 

In this paper, we study Hamiltonian lattice QCD with Wilson fermions at finite chemical 
potential. We derive the effective Hamiltonian in the strong coupling regime and diagonalize 
it by Bogoliubov transformation. The vacuum energy, chiral condensate, and masses of 
pseudo-scalar, vector meson and nucleons are computed. In the non-perturbatively defined 
chiral limit, we obtain reasonable results for the critical point and some physical quantities 
in the large N c limit, with N c the number of colors. 

To our knowledge, the only existing literature about the same system (r^O and /i ^ 0) is 
Ref. [21], where the author used a very different approach: the solution to the gap equation. 
In contrary to the conventional predictions [19], the author found that even at fx — 0, there 
is a critical value for the effective four fermion coupling K, below which dynamical mass 
of quark vanishes. He introduced the concept of total chemical potential and found the 
transition order depends on the input parameters K and r as well as the momentum. In 
contrast, we find that in the chiral limit, dynamical mass of quark doesn't vanish for all 
values of K if ji < ji c (the chiral-symmetry broken phase); and at /i — /ic, our order of 
chiral phase transition doesn't depend on the input parameter. 

The rest of the paper is organized as follows. In Sec. II, we derive the effective Hamiltonian 
at finite chemical potential. In Sec. Ill, we present the results for the vacuum energy, chiral 
condensate, and hadron masses. In Sec. IV, we estimate the critical chemical potential at 
zero temperature. The results are summarized in Sec.V. 
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II. EFFECTIVE HAMILTONIAN IN THE STRONG COUPLING REGIME 

A. The /j, = case 

We begin with QCD Hamiltonian 1 with Wilson fermions at chemical potential fx = on 
1 dimensional continuum time and 3 dimensional spatial discretized lattice, 

1 _ ±d 

H = M Y,$(x)if>(x) + - E E $(?)lkU{x, k),/>(x + k) 

x Za x k=±l 

±d _ 2d 

" orE E j(x)U(x,kMx + k) + ^'£E?(x)E?(x) 

Za x k=±l Za x j=l 

-iETr(^ + f/ p + -2), (2.1) 
u y P 

where 

TO 

M = m + —, (2.2) 
a 

d — 3 is the spatial dimension and m, a, r and g are respectively the bare fermion mass, 
spatial lattice spacing, Wilson parameter, and bare coupling constant. U(x, k) is the gauge 
link variable at site x and direction k. Fermion field tp on each lattice site carries spin 
(Dirac), color and flavor indexes; here and in the following, whenever there is a summation 
sign "X)x") summations over these indexes are implied. The convention = — 7^ is used. 
E"(x) is the color-electric field at site x and direction j, and summation over a = 1, 2, 8 
is implied. U p is the product of gauge link variables around an elementary spatial plaquette, 
and it represents the color magnetic interactions. In the continuum limit a — > 0, Eq.(2.1) 
approaches to the continuum QCD Hamiltonian in the temporal gauge A 4 = 0. 

The effective Hamiltonian, obtained by strong coupling expansion up to the second order, 
is [20] 



We use a representation of 7 matrices described in D. Lurie, Particles and Fields, Interscience 
Publishers, John Wiley & Sons, Inc. (1968). 
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+ (r 2 - l)^ h (x)^ h {x)^) 2 {x + k)^4) h {x + k) 

- (r 2 - 1)^(^)75^/2 (x)i/j\ 2 (x + k)^ fl (x + k) 

+ (r 2 + l)^} i (x)7475V' /2 (x)^} 2 (x + k)lilh^hi x + fc) 

+ (r 2 + (1 - 25| fc |j)) ^(^747^/2 WVi^ + k)l±l^ h (x + fc) 

- (r 2 - (1 - 2<J| fc |j)) ^(^7^/2(^)^/2 (a: + k)^ fl (x + k) 

- (r 2 + (1 - 25\ k \ d )) ip^ix^Vjipfrix^frix + k)-f 4 aji) fl (x + k) 

~ (r 2 - (1 - 2S WJ )) ^ fl (x)a^ h (x)^ f2 (x + k)a^ h {x + fc)). (2.3) 



This effective Hamiltonian is equivalent to that in Ref. [19], though different representations 
of 7 matrices are used. Here Oj = £jj 1 j 2 ljilj2- The flavor indexes are explicitly written. This 
effective Hamiltonian describes the nearest-neighbor four fermion interactions, with 



being the effective four fermion coupling constant. Here Cjv = (iV 2 — 1)/ (2N C ) is the Casimir 
invariant of the SU(iV c ) gauge group. 

B. The fj, ^ case 

In the continuum, the grand canonical partition function of QCD at finite temperature 
T and chemical potential ji is 



Z = Tre-M H -' tN) , (3={k B T) 



-1 



(2.5) 



where k B is the Boltzmann constant and N is particle number operator 



AT = Y,^(x)ij(x). 



(2.6) 



X 
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According to Eq.(2.5) and following the procedure in Sec. II A, the role of the Hamiltonian 
at strong coupling is now played by 

Heff = H eff - ptN. (2.7) 

In this Hamiltonian, there are three input parameters: r, m and ji. Suppose we study the 
phase structure of the system in the chiral limit. Such a limit can be reached by fine-tuning 
the bare quark mass m so that the pion becomes massless. In such a case, there are only 
two free parameters left: r and //. 

The vacuum energy is the expectation value of H — /iN in its ground state and also 
the expectation value of H^fj m ^ s g roun d state |^ e //), given by 

E n = (Sl\H- pN\Q) = (n eff \H» ff \n eff ). (2.8) 



III. PHYSICAL QUANTITIES AT fx / AND T = 
A. Meson masses 

One way to compute the masses of mesons as well as the contributions of the mesons to 
the vacuum energy, is to bosonize the effective Hamiltonian Eq. (2.3). We introduce the 
following operators [9,22,23] 



n /l/ 2 (^) = 


2y 


1 


^\{x){\ 74)75^/2(3), 

-V 


11/2/1 ( X ) = 


2, 


1 


tf h {x){l + 74)75^/1(3), 


YifiM = 


2, 


1 


—^ h (x)(l 74)7^/2(3), 




2a 


1 


-tp} 2 (x){l + 74)7^/1(3). 



j stands for the positive spatial direction, and v and v* denote respectively expectation value 
of ipip and ip^ip in the vacuum state |^ e //) of H e ff, i.e., 
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I eff ^f^s 
1 



\ /eff N f N s 

Here N s is the total number of lattice sites and Nj the number of flavors. It is shown in 
Appendix A that under the linearization prescription [23], operators defined in Eq. (3.1) 
satisfy the canonical commutation relations for bosons and then the effective Hamiltonian 
H£jf in Eq.(2.7) can be expressed in terms of these operators in the following way 



where 



Hfa ~ m mear = 4 0) + H u + H V , (3.3) 

„(o) » r »J„ (Kd(l + r 2 ) \ t Kdr 2 , + x 
4 0) = N f N s \Mv - f K - + // W + — ( V J + u 2 ) 



and 

\ c / X,/l,/2 

c x,fi,h,k 
\ c / x,f 1: f 2 ,j 

Kr 2 

+ ^ E K/ 2 (*) (* + k) + V jhh {x)V} hh (x + kj) 

c x,fi,f 2 ,k,j 

~ 4aA^ S (^U^)4 2 / 1 ^ + ^) + % 1 / 2 ^)% 2 / 1 ^ + ^))( 1 - 2 ^)- ( 3 - 5 ) 

c x,fi,f 2 ,k,j 

vl, V2, v\ c . and are expectation value of four fermion operators 
^2 = * AT ^ e// * ^ i ) hi x )^f2{x)^ h {x + (x + k)\n eff ), 

f s x,fi,f 2 ,k 
V2 = 9JM AT ( Qe ff\ E ^(^V^W/^ + k)^ fl (x + k)\Q eff ), 

zaiyiM s x ,h,f 2 ,k 



f S X,fl,f2,k,j 



2dN f N 8 x ,-r fc ■ 



Hn in Eq.(3.5) becomes 



aiV, 



a-*- 2 )) E n} l/2 (p)n /2/l (p) 

> p,h,h 



The Bogoliubov transformation [20] 



n(p) — > n(p) cosh -Up + Ir(— p) sinh 



(3.6) 



After a Fourier transformation 

n /l/2 (x) = £e^n /l/2 (p), (3.7) 



+ oTat^ E E (n} l/2 (p)n /2/l (p) + n(p) /l/2 nj 2/i (p)) £ C os Pj a 

Z " iVc /l/2 V j=l 

~^N d E K A (p)nU(-p) + nAA(-p)nAA(p))E«»Pia- ( 3 - 8 ) 

c p,/i,/2 i=i 



IF(p) — > IF(p) cosh -Up + n(— p) sinhM p , (3.9) 

diagonalizes ifn if 

tann 2-Up = — — — 7 y cosp/a, 
G i j=i 

Gi = 2M - -^(1 - r 2 ) + — — u ^ cosp,a, 
aN c aN c fr[ 

G 2 = —v. (3.10) 

2aN c K ' 

The resulting H u is 

H U = G 1 E (1 - tanh 2 2 M p)^n} i/2 (p)n /2/l (p) 

P,/l,/2 

" ^/E ((l " (1 " tanh 2 2^) + H^£cosp,a) . (3.11) 



Now n^y 2 (p) stands for the pseudo-scalar creation operator in the momentum space. Ac- 
cording to Eq. (3.11), the difference between the pseudo-scalar meson energy and vacuum 
energy is 



E n = Gi(l- tanh 2 2u p ) 



= [iM - - r 2 ) - -^-(1 - r 2 )v £ cos Pl a^ 



i 

d \ 2 



(Kd K 
2M- ^iV t;(1_r2) + ^V ( 1 + r2 )^J2 C0S Pi a ) > ( 3 - 12 ) 

which gives the pseudo-scalar mass when pj = 0. The pseudo-scalar mass square is 

In order to define the chiral limit, one has to fine tune M — > M chira i so that the pion becomes 
massless. From Eq. (3.13), we get 

Kdr 2 

M chiral = —v. (3.14) 

aN c 

In this limit, the pseudo-scalar mass square behaves as oc M — M c ^ m /, which is the 
PCAC relation. 

The Hy sector in Eq. (3.5) can be considered in a similar way. After a Fourier transfor- 
mation 

V 3 {x) =^%{p) (3.15) 
v 

and a Bogoliubov transformation 

Vj(p) -> Vj(p) coshw { p j) + V-(-p) sinhw^, 

Vj(p) — > t^'(p) coshwj(p) + Vj(—p) smhiu^, (3.16) 

H v becomes a diagonalized one 

H V = G 1 £ (l-tanh 2 2^))^^ l/2 (p)% 2/l (p) 
P,3,h,h 

- ^/E ((l " (l - tanh 2 2^)) 1 ) + ^r 2 fcos^ , (3.17) 



if 



2G2 d 

tanh 2-u^ = — (E cos ^ a — 2cospjd). (3.18) 
Gi 



1=1 

T/t 



Now Vjfafeip) stands for the vector creation operator in the momentum space. The vector 
mass is 



M v = G 1 (l - tanh 2 2^ } ) 1 — ^m—*m c _^X^Jy. (3.19) 
According to Eqs. (3.3), (3.11) and (3.17), the vacuum energy reads 

E n = (n\H»\n) 

= - ^iV 2 E ((l - (1 - tanh 2 2urf) + ^ E cos W a) 

- ^E (f 1 " (1 - tanh 2 2^'))^) + ^r 2 ]T cos ^ . (3.20) 

1 P,3 \ Gl 1=1 / 

This also gives the thermodynamic potential (grand potential) at T = 0. 

As shown in Ref. [20], at /i — the results above are consistent with those of Smit in 
Ref. [19] where the 1/N C expansion was used. 

B. Results in the large N c limit 

As shown in Appendix B, in the chiral limit, the dominant contributions to the vacuum 
energy for large N c is 

~2\ \ ZTJjl 



E Q -> NfN s 



(Kd(l + r 2 ) \ , Kdr 2 ( , A 2 „\ Kd ( 2 / A 2\ 



(3.21) 



As shown Appendix C, under the mean-field approximation, i.e., by Wick-contracting 
a pair of fermion fields in the four fermion terms in Eq. (2.3), one can obtain a bilinear 
Hamiltonian in the large N c limit 

H eff ~ H MFA = A^W(x) + BY,i>\x)1>{x) + C, (3.22) 
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where 

A = M chiral - ——{1 - r 2 )v, 
2al\ c 

B = a {Wr 1 )' 

C = ((1 + r V 2 - (1 - r V) N s N f . (3.23) 

The coefficient A plays the role of dynamical mass of quark. According to Eq. (2.8) and Eq. 
(3.22), the vacuum energy in presence of the chemical potential \i is now Eq = (Q e ff\H MF A — 
/j,N\Q e ff), which agrees with Eq. (3.21), derived from the large N c limit of a bosonized 
Hamiltonian Eq. (3.3). 

In the large N c limit, the fermion field ip can be expressed as 



(3.24) 



The 2-spinors £ and Ty are the annihilation operator of positive energy fermion and creation 
operator of negative energy fermion respectively. Let us define the state \n p ,n p ) in the 
momentum space by 

£p|0p, Up) — 0, ^p|0p,77,p) = |lp, Tip), ^p\lp,flp) = I Op, Tip), ^p|lp,7T,p) = 0, 

T] p \n p , Op) = 0, vt\ n p, Op) = \n p , l p ), 7] p \n p , l p ) = \n p , Op), ^|n p , l p ) = 0. (3.25) 

The numbers n p and n p take the values or 1 due to the Pauli principle. By definition, the 
up and down components of the fermion field are decoupled in such a state \n p , n p ). For the 
vacuum state of H^ffi we m& ke an ansatz 

Peff)= fn p ,n p \n p ,n p ). (3.26) 

In the leading N c limit, the dominant contributions to the chiral condensate and quark 
number density are 

WV/ N f N s 
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S C npfip (n p ,n p \^\n p ,n p ) = 2N c (n + n-l), 

n p ,n p ,p 



N f N s 



n q 



{n\Y. x ^{x)j;{x)\n) 

2N c NfN s ~* 2N C 

= 0AT AT AT S C n P ,n p (n p ,n p \ij^\n p ,n p ) -l = n-n. (3.27) 

Here we denote C„ Pi n p = fn p ,n p - Using Eq. (3.21) and Eq. (3.27), in the large N c limit we 
obtain the normalized vacuum energy 

en ~ 2N c N f N s 

Kdr 2 

= M chira i (n + n-1) (n-n + 1) 

a 

Kdr 2 too \ Kd , . . 

H [n + n + 1 — 2n) -\ (n + n- 2nn - 1) - fi (n - n + 1) , (3.28) 

a v 'a 

where the quark number n and ant i- quark number n 

n= (n) = C n P ,n p n p , 

n p ,n p ,p 

n= (n) = C n P ,n p n P (3.29) 

n p ,n p ,p 

are constrained in the range of [0, 1] and determined by minimizing the vacuum energy. 
For a generic nucleon operator O^uci consisting of three quarks, the thermo mass is 

M Nud = (VL efS \0 Nud H^ ff O ] Nud \VL ef f) - E n . (3.30) 

Under the mean-field approximation in the large N c limit, it becomes (see Appendix D) 

M Nucl = 3{A + B) - 3/i. (3.31) 



IV. PHASE STRUCTURE AT T = AND fj, ^ 

We now consider in the larger N c limit, the critical behavior of the system at T = 
and /i ^ 0. The ground state of the system corresponds to the lowest value of the vacuum 
energy. At some given inputs of Wilson parameter r and chemical potential /i, we can find 
the value of n and n when Eq. (3.28) is minimized. The result is 

12 



n = 0(/i - nc), 

n = 0. (4.1) 

Here G(/x — He) is the step function: it is for /j, < jj, c and 1 for /i > jj, C) where He is the 
critical chemical potential 

l*c = — (l + 2r 2 ). (4.2) 

(X 

Substituting Eq. (4.1) into Eq. (3.27), we obtain the chiral condensate and quark number 
density 

{W) = 2N c (e ^-^ c )-l), 

n q = e(fi-fx c ). (4.3) 

There is clearly is a first order chiral phase transition. For \i < /ic, the system is in the 
confinement phase with chiral-symmetry breaking. For /i > fi C , chiral symmetry is restored. 
According to Eqs. (3.23), (3.31) and (4.1), the thermo mass of the nucleon is 

M Nud = M$l cl - 3/1, (4.4) 

where 

MSL = 3m£ (4.5) 

is the nucleon mass at /i — 0, and 

m dyn = ( 4 -6) 

is the dynamical mass of quark at /j, — 0. From Eq. (4.4), one sees that the nucleon thermo 
mass vanishes at [i — M^ ud /3 = mfy n , before the chiral phase transition takes place. This 
is not surprising because Wilson fermions break explicitly the chiral symmetry. The value 
of fie should coincide with M^ ucl /3 when r is very small, i.e., as in the case of naive or 
Kogut-Susskind fermions [9]. 
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V. DISCUSSIONS 



In the preceding sections, we have investigated (d+l)-dimensional Hamiltonian lattice 
QCD at finite density with Wilson fermions in the strong coupling regime. We compute the 
vacuum energy, meson and nucleon masses, chiral condensate and quark number density. 
At finite chemical potential, there is an interplay between the bare fermion mass in the 
chiral limit and the chiral condensate, which has to be determined self-consistently. The 
critical behavior of the system in the large N c limit is considered: a first order chiral phase 
transition is found at ji = fie] the nucleon thermo mass vanishes before the chiral transition 
takes place, which is due to the explicit breakdown of chiral symmetry by Wilson fermions. 

We have not yet specified the nature of the chiral- symmetric phase for fj, > /ic- Is it a 
QGP phase or a color-superconducting phase [24,25]? Up to now, there has been no first 
principle investigation of such a phase in SU(3) gauge theory. The answer to this question 
might be very important to our understanding of the formation of the neutron star or quark 
star. 

We also know that the strong coupling regime is far from the continuum limit. One 
has to develop a new numerical method to study the continuum physics. The Monte Carlo 
Hamiltonian method developed recently [26,27] might eventually be useful for such a purpose. 
We hope to discuss these interesting issues in the future. 

ACKNOWLEDGMENTS 

We thank V. Azcoiti, S.H. Guo, H. Kroger, V. Laliena, and M. Lombardo for useful 
discussions. X.Q.L. is supported by the National Science Fund for Distinguished Young 
Scholars (19825117), Key Project of National Science Foundation (10235040), Guangdong 
Natural Science Foundation (990212), National and Guangdong Ministries of Education, and 
Hong Kong Foundation of the Zhongshan University Advanced Research Center. Y.Z.F. is 
partly supported by the National Natural Science Foundation (10235040 and 10275098). 



14 



APPENDIX A: BOSONIZATION AND LINEARIZATION 



In Sect. Ill, to bosonize the Hamiltonian Eq. (2.3), we introduce the operators II, 11*, 
V and V*, defined by Eq. (3.1). In Eq. (2.3), there are terms having direct correspondence 
to these operators. In addition, there are also terms irrelevant to these operators. 

In order that the operators in Eq. (3.1) represent the appropriate mesons, they must 
satisfy the commutation relations for boson operators. However, a direct calculation shows 



n /i,/ 2 (^),n} 2j/l (x)] = — {4) h {x)i) h {x)+i) h {x)il) h {x)) , 



-j= (^h{ x )(liH ~ lHj)M x ) + iPfr{ x )(ljli - Hlj)^h{ x )) , ( A1 ) 



2v 

2v 
1 

Av 

which are not consistent with the commutation relations between the annihilation and cre- 
ation operators for bosons. Similar situation also appears in quantum theory of magne- 
tization [28] and many-particle systems [29], where a linearization prescription is used to 
simplify the theory. Using such a procedure, i.e., with the fermion bilinears in the r.h.s. of 
the commutation relations replaced by the vacuum expectation value ip(x)ijj(x) — > v and for 
degenerate quarks, Eq. (Al) becomes 

[U(x),U (x)\ ~ - — - 1, 

[Vj{x), V, [x)\ ~ d fl j^- - d fl , (A2) 

which are the correct commutation relations for bosons. As shown in Refs. [20,23], the 
linearization prescription leads to results consistent with Ref. [19], where systematic 1/N C 
expansion was used. 

The four fermion terms in Eq. (2.3), which have direct correspondence to the boson 
operators, are 

H Wi ( x h^f 2 (a#/ 2 ( x + fc)75^/i ( x + k) + V'/i {x)iii^f2 ( x )^f 2 ( x + fc)7475^/i (x + k)) 

= - 2 ^E K/>) n L> + *) + n /l/2 (x)n /2/l (x + k)) , 
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E (^/i (a075^/ 2 (a#/ 2 (x + fc)75^/i (x + fc) — ^ (a;)7475^/ 2 (x)4>\ 2 (x + k)^^ h (x + fc)) 

= -2vZ K/»n /2/l Or + fc) + n /l/2 (x)nt 2/i (x + k)) , 

x,k 

E (^/i ( x )747i^/ 2 (a#/ 2 (a: + k^lj^h (x + k) + ^ (a:)7j^/ 2 (a#/ 2 (x + khjiph (x + fc)) 

x,k,j 

x (l - 25| fc |j) 

= -2* E fcW^W* + *) + + *)) (l - , 

E Wi Wni^h (a#/ 2 (a: + ~k)lilj^h (x + k) - ^ (x)7^/ 2 (x)^} 2 (x + fc)7j^/i (s + fc)) 

(A3) 



x,k,j 



^ E {v} hh (x)V jhfl (x + k) + V jflh (x)V} f2fl (x + k)) . 

x,k,j 

For the bilinear J2 X ^{x)%l){x) in Eq. (2.3), a direct calculation shows 



1 



J2^(x')ip(x'),U flh (x) 



V'} 1 (x)7 4 (l -74)75^/2^) 



-^(^74(1 - 74)7574^/ 2 (a;) 

^=^} l (x)(i - 74)75^/2^) = -n /l/a (x). 



Similarly, one can also show 

"1 



2J-v 



J2i>(x'Mx>),n\ 2fi (x) 



(A4) 



-^(x'Mx'),V jhh (x) 



= ~ V jflf2(x), 



^E^W),4 2 />) =4aW. ( A5 ) 

x' 

Eq. (A4) and Eq. (A5) imply that the bilinear ^2 x ip(x)ip(x) can be bosonized as 

J2^Mx) ~ E f^ + 24 2/i (x)n /l/2 (x) +2^4 2/l (x)% l/2 (x) J . (A6) 
For the four fermion term Y,x'4 ) h{ x )' l l ) h{ x )'4 ) h{ x + ^)^fi( x + ^) m Eq. (2.3), we have 

\ t n (x')^ W + W + k) , n /i/ 2 (*) 

x' 

E ^ W/ 2 (*') \l$r 2 (x' + fc)^ (*' + k) , n /l/2 (x) 

x' 

■1 — 



+ E ^/^W*')' 11 /!/^) $ f2 (x' + k)^ n (x' + k) 

x' 

- -2^ n /i/ 2 (^), 



(A7) 
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where in the r.h.s., Eq. (A4) and the linearization prescription have been used. In addition, 
we also have 



\ £ ^ (x')^ {x')$ f , (x' + k)^ (x' + k) , 4 2/i (x) 

x' 

- J2^f[(x')^(x')^(x' + k)ij f[ (x' + k), V jflh (x) 



1 



- H V)$f 2 ( x> + k)iP fl (x' + k), V} hh . 



X) 



- -1vV jhh {x), 
~2vV} Ml {x). 



(A8) 



2 

x' 

Eq. (A7) and Eq. (A8) also imply that the four fermion term J2x4>fi( x )i J f2( x )i'f2( x + 
k^f^x + k) can be bosonized as 



E i>h ( x + k)4>A ( x + k) 

X 

£ i) 2 + 44 /i (x)n /l/2 (x)+4^4 /i (x)% /2 (x) ) . 



(A9) 



In the linear prescription, the four fermion operators in Eq. (2.3), which are irrelevant 
to the operators in Eq. (3.1) are replaced by their vacuum expectation values. For example, 

ip\x)ip(x) ~ v\ 
E A ( x )^h i x )A ( x + kWfi ( x + k)~ 2dN f N s v\, 
E (r 2 + (l - 2<S| fc |j)) $fA x )o$h( x )$h( x + tyv^hix + k)~ 2dN f N s v 2 a j , 



x,k,j 



E (? - ^ - A( x ) a ^h( x ^h( x + k)a^ h {x + k)~ 2dN f N s vl.. 

x,k,j 

Collecting the above results, we obtain Eqs. (3.3), (3.4) and (3.5). 



(A10) 



APPENDIX B: VACUUM ENERGY IN THE LARGE N c LIMIT 

According to Eq. (3.6), 

^ 2 = 9^f AT ^ ((^fuM^h,^)) (i ) f2,c 2 { X + k)*pf UC2 ( x + k)\ 

2dN f N s x kJij2 V\ /eff\ I eff 



+ \ i , hi. x )i , f2,c 2 { x + k) ) (ip h (x + k)ijj fuci (x) 

eff \ / eff 
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2 



+ + ty eff \M x + ^^.- (x) / e J' (B1) 

where the color and flavor indexes for the fermion fields are explicitly specified and summa- 
tion over the repeated color index is implied, (--^eff stands for the expectation value taken 
in vacuum state \Q e ff). Because 

ijfi,ci(x)^f 2 , ci (x) 





= v5 h,h 


eff 






<5ci,c 2 i 


eff 






= V ^h,h 


eff 






^C1,C2 ? 


eff 





1 ' - - m \ ■ - 'B2) 



terms from Wick-contraction of bilinear at different lattice sites are subdominant for large 
N c . Therefore, we have 

V2~^v 2 , 4 — (B3) 

with v and expectation value of fermion bilinears defined in Eq. (3.2). Expectation values 
of other four fermion operators in Eq. (3.6) can be considered in the similar way, 

v 2aj oc (Tr^) 2 = 0, vl. oc (Traj) 2 = (B4) 

in the large N c limit. 

In the chiral limit Eq. (3.14), G 1 defined in Eq. (3.10) scales as vK/N c oc K and G 1 /G 2 
doesn't depend on N c . The Nj terms in Eq. (3.20) are subdominant because they are 
proportional to K, while other non- vanishing terms of E$ in Eq. (3.20), defined in Eq. 
(3.4) are dominant because they behave as vK oc KN C or v^K oc KN C . Using this fact and 
substituting Eqs. (B3) and (B4) into Eq. (3.4), we obtain Eq. (3.21) from Eq. (3.20), i.e., 
the vacuum energy in the large N c limit. 
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APPENDIX C: MEAN-FIELD APPROXIMATION IN THE LARGE N c LIMIT 



The mean-field approximation is a popular technique wildly used in quantum field theory 
with four fermion interactions (for example, to bilinearize the Nambu-Jona-Lasino model 
[30,31]) and quantum theory of many-particle systems [29]. 

In order to show how to derive Eq. (3.22) from Eq. (2.3), let's look at the first and last of 
the four fermion terms in Eq. (2.3) as an example. Under the mean-field approximation, i.e., 
replacing a pair of fermion fields by their vacuum expectation value, the first four fermion 
term in Eq. (2.3) becomes 

x,k,fi,f 2 

~ E (Muci( x )^h,ci( x )) ^} 2 ,c 2 ( x + k)ip fl ,c 2 {x + k) 
x,k,h,h v x 1 eff 

+i>\ ucl (x)4>h,cA x )(i>} 2 ,c 2 ( x + k)ipf UC2 ( x + ty) 

+^h,ci( x )(^hA x )^f 2 ,c 2 ( x + k)) ip fl , C2 {x + k) 

X ' eff 

+^h,c 1 (x)^ f2 , C2 ( x + k)(^} 1:Cl (x)^ fl , C2 (x + fc)^ ) - 2dN f N s vt (CI) 

Again, the color and flavor indexes for the fermion fields are explicitly specified and sum- 
mation over the repeated color index is implied. According to Eq. (B2), for large N c , Wick 
contractions of bilinear at different lattice sites give subdominant contribution, therefore 

E i'f 1 ,c 1 ( x )^f2,c 1 ( x )^f 2 ,c 2 ( x + k)ip fl , C2 ( x + k) 

x,k,fi,f 2 

~ E 5 h ,h (V} 2 > + k)ip h , c ( x + - 2dN f N s (v^f 

x,k,h,f 2 V J 

= « f E (i>U x )M*) + i>U x + + *)) - 2dN f N ° ( vt ) 2 

x,k,f 

= - 2dN f N s (V) 2 . (C2) 

X 

Under the mean-field approximation in the large N c limit, the dominant contribution to 
the last four fermion term in Eq. (2.3) is 
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~ H [(^fucii^j^h,^)) ^h,c 2 ( x + k)<Tjil;f u0a (x + k) 

x,k,fi,f 2 V X 1 e ^ 

+^f 1 ,c 1 (xyjiPh,c 1 (x)(ip} 2!C2 (x + k)<Tjipf llCa (x + k)) ) - 2dN f N s vl a . 

x ' e/// 

OC 7V(<7j) = 0. (C3) 

Treating other four fermion terms in the same way, one obtains Eq. (3.22), i.e., the bilinear 
Hamiltonian Hmfa- 



APPENDIX D: NUCLEON MASS 

Under the mean-field approximation in the large N c limit, Eq. (3.30) becomes 

M Nucl = (n eff \0 Nucl (H MFA - fiN) ] Nucl \tt eff ) - E n . (Dl) 

Substituting the bilinear Hamiltonian Eq. (3.22), the fermion field Eq. (3.24), and the 
particle number operator Eq. (2.6) into Eq. (Dl), we have 

M Nucl = (n eff \0 Nucl Ia^H x )H x ) + (B- ii)Y t ^\x)rl){x) + c\ 0^ Nucl \Q eff ) - E n 

\ X X / 

\ X X / 

= (n eff \0 Nud (AJ2et + (B-»)j:^) O j Nuc iKff), (D2) 

\ X X / 

where we have used the fact that the terms with rfrj and C are canceled by En- 

Let us take the proton as an example. One can write the operator 0^ Nucl explicitly as 

O^Nud = ^gjy X! e c 1 c 2 c 3 ^l 1 ,u,l( X ) (^,l 2 ,u,l( X )^l 3 ,d,2( X ) ~~ £,i 2 ,u,2( X )£,i a ,d,l(. X )) j (D3) 

where Ci,c 2 and c 3 are the color indexes, u and d stand respectively for the u-quark and 
d-quark, and 1 and 2 are the spin up and down indexes. Using the anti-commutation 
relationship for fermions 

ic,f, s {x)il l fl sl {x') = 5 c>cl 5 fJ ,5 S:S ,5 x/x , - Cf lt Ax')^ fta (x), (D4) 
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one obtains 

E^WOkA//) = SOkdl^e//). (D5) 

x 

Eq. (3.31) is a consequence of Eq. (D2) and Eq. (D5). 

At 11 — 0, the nucleon mass equals to 3m^ n (see Eq. (4.5)), which agrees with Refs. 
[19,32]. 
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